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We study the properties of cosmological density perturbations in a multi-component system con- 
sisting of a scalar field and a perfect fluid. We discuss the number of degrees of freedom completely 
describing the system, introduce a full set of dynamical gauge-invariant equations in terms of the 
curvature and entropy perturbations, and display an efficient formulation of these equations as a 
first-order system linked by a fairly sparse matrix. Our formalism includes spatial gradients, ex- 
tending previous formulations restricted to the large-scale limit, and fully accounts for the evolution 
of an isocurvature mode intrinsic to the scalar field. We then address the issue of the adiabatic 
condition, in particular demonstrating its preservation on large scales. Finally, we apply our for- 
malism to the quintessence scenario and clearly underline the importance of initial conditions when 
considering late-time perturbations. In particular, we show that entropy perturbations can still be 
present when the quintessence field energy density becomes non-neghgible. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 



The material content of the universe is commonly as- 
sumed to be a mixture of fluids, such as radiation or 
non-relativistic matter, and scalar fields, either driving 
a period of early universe inflation [IJ or playing the 
role of dark energy (quintessence) in the present universe 

1^ Q . The latter possibility has motivated a number 
of works devoted to the study of cosmological pertur- 
bations in a multi-component system consisting of fluids 
and a scalar field, for instance Refs. 

0,0. Nevertheless, the literature contains some con- 
tradictory statements concerning the properties of such 
perturbations. 

In this paper we aim to resolve these discrepancies and 
will provide a comprehensive analysis of the problem. We 
will study the role of intrinsic entropy perturbation in the 
scalar field, and whether the notion of adiabaticity is pre- 
served by the dynamics of the multi-component system 
when the evolution of such an intrinsic entropy pertur- 
bation is explicitly accounted for. In the process of do- 
ing so, we will discuss the number of degrees of freedom 
which completely describe the system and we will find a 
highly-efficient formulation of the perturbation equations 
including the effects of spatial gradients. 

Finally, using our formalism we will specifically discuss 
the quintessence scenario. We will correct some common 
misconceptions, discuss the evolution of entropy pertur- 
bations, and clearly show the importance of initial condi- 
tions when considering late-time perturbations. In par- 
ticular, we will show that entropy perturbations can be 
enhanced by the evolution of the field and may still be 
present when its density is no longer negligible. This is 
an important result which has generally been overlooked 
when studying dark energy models. 
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II. THE DYNAMICAL EQUATIONS 

A. Background 

Our approach builds on an earlier paper by Malquarti 
and Liddle jjl, and we will largely follow the notation 
of that article but with some differences in definitions. 
We assume a flat Friedman-Robertson- Walker universe 
throughout, with the background evolution determined 
by the usual equations 



2H + 3H^ = 2- 

a 



STrGptot , (1) 
-S^Gptot , (2) 



where H = d/a is the Hubble parameter, a is the scale 
factor, a dot stands for a derivative with respect to cos- 
mic time t and the subscript "tot" always refers to the 
sum over all matter components. The fundamental ingre- 
dients we consider here are a perfect fluid with constant 
equation of state Wf = Pf/pf and a minimally coupled 
scalar field with potential V{ip). Since we treat the 
fluid and the scalar fleld as uncoupled, the conservation 
of their respective energy-momentum tensors gives 



Pf 



dip 



-3(l+«;f) 



(3) 
(4) 



The subscripts "f" and "1^9" will always refer to the per- 
fect fluid and the scalar fleld respectively. 

Useful parameters describing completely the back- 
ground properties of the scalar field are its equation of 
state 



Pip 
Pv 



ip'^/2 + V{ip) 



(5) 



2 



and its adiabatic sound speed 



w, 



2 _ 'Pv_ _ _ 



= 1 



2 dV/dip 

3 Hip 



(6) 



We also introduce the total equation of state w = 
Ptot/ptot and the total sound speed Cg = ptot/ptot, and 
in order to simplify some expressions we define for each 
component 7x = 1 + and also 7 = 1 + 



B. Perturbations 

We consider only scalar perturbations and we choose 
to work in Newtonian gauge where the perturbed 
metric reads as 



ds' 



(1 + 2$)di2 + a'^{t){l - 2*)dx2 



(7) 



Here 4> and 4" describe the metric perturbation, and in 
this case are equal to the gauge-invariant potentials de- 
fined in Ref. |l3j | . We work in Fourier space and compute 
the first-order perturbed Einstein equations. As our sys- 
tem has no anisotropic stress, the (i — j) Einstein equa- 
tions imply that the metric potentials are equal, 4* = $. 
The remaining Einstein equations are 



3i/(i/$ + $) - 



{2H + 3H^)^ 



47rG<5ptot , (8) 

AirGSptot , (9) 
47rG<5<Ztot , (10) 



where VJgtot is the total momentum perturbation of the 
system. Eq. © comes from the (0 — 0) Einstein equation, 
Eq. (jSJ from the (i—i) Einstein equation, while Eq. (|10|1 is 
obtained from the (0 — i) Einstein equation. The perfect 
fluid and scalar field perturbation variables are 



6pf 


= WfSpf , 




(11) 


5qi 


= Pf7fVf , 


dp 


(12) 




= ipSip — p^^ + 


(13) 




= ipSip — — 


dp 


(14) 




— —ip5p, 


(15) 



where Vt is the fluid velocity potential defined so that 
the fluid velocity is given by 5vl = VVf — this is pos- 
sible since for scalar perturbations the flow is irrota- 
tional. Note that this definition is slightly different 
from the one used in Ref. The conservation of the 

energy-momentum tensors for each component provides 
the equations 

1,2 A2y Ay 

Sip + 3HSp + —Sp + -^5p = 4(^$ - 2—$ , (16) 



dp^ 



dp 



5i — 37f<l> 



7f^ 



Vf, 



Vf = mwfVf — -5f-^ 

7f 



(17) 
(18) 



where 5f = 5pi/pi. We also define 5^ = Sp^jp^p and 
5^ = ^v^pl p^p- These equations are the perturbed Euler- 
Lagrange equation for the scalar field and the continuity 
and Euler equations for the fluid. 

As we will see later, it is useful to introduce the co- 
moving density perturbation for each component |l4j 



(19) 



which is a gauge-invariant quantity. We also introduce 
the gauge-invariant entropy perturbation variables [T^L 
Il6lll7l |. namely the relative entropy perturbation between 
the fluid and the scalar field 



S . 



3-ff7f7ipilf ( 5p^ Spi 
7 V Pv Pi 



lA^Jl^, (20) 

7 



and the intrinsic entropy perturbation of the scalar field 

5-K 



1 



Pv 



Pv 



1 - c2 



(21) 



where Stt = 5p^p/ p^p. The normalizations have been cho- 
sen in order to simplify some later expressions. Note that 
r is well defined even if c^^ ~ 1 since, as can easily be 
shown, we have T = t^j p^. By definition the perfect fluid 
does not have an intrinsic entropy perturbation. Adia- 
baticity is defined by the condition S* = F = 0, since in 
this case it is possible to define a slicing for which all 
matter component perturbations vanish. 



C. Degrees of freedom 

The system of Eqs. ©, ®, 119, (113, (113 and ^ 

describes the evolution of 4 variables, namely $, pf, Vf 
and bp. Equations © and (|16|l are second order and if we 
introduce two new variables for $ and bp (and therefore 
two new equations) we end up with 6 variables describ- 
ing the perturbations, 6 first-order dynamical equations 
and two constraint equations (Eqs. ((HJ and HlOfl ^. The 
two constraint equations reduce the number of degrees 
of freedom to 4 and as a result two dynamical equations 
must be redundant. This comes from the fact that the 
conservation of the total energy-momentum tensor is a 
consequence of the Einstein equations and therefore the 
conservation equations for one matter component implies 
the ones for the other. As a result, it is possible to 
write this system as four differential equations for the 
four dynamical degrees of freedom which completely de- 
scribe the perturbations, e.g. <&, bi^^^o and b-,^ which are 
the basic variables studied in Ref. |ll| in the large-scale 
limit k/aH ^1. A general solution to those equations 
permits both an isocurvature perturbation between the 
scalar field density contrast and the fiuid density con- 
trast, and an isocurvature perturbation intrinsic to the 
scalar field, i.e. between its density and pressure pertur- 
bations. 
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For our discussion, it is useful to combine Eqs. © and 
(|10|1 and find the constraint equation 



fc2 

$ 

a2 



-47rGeb 



(22) 



Note that this is a gauge-invariant equation, though had 
we included an anisotropic stress then $ must be replaced 
by the second metric potential ^' 0|. If the fluid is 
completely absent, so that we simply have a single scalar 
field, the constraint equation Eq. (|22|l reduces to 



-$ = 



47rG 
172" 



(p5^)^--T. (23) 



The system is completely described by two dynamical 
degrees of freedom and this equation implies that one of 
the scalar field degrees of freedom is removed, e.g. Sip. 
The right-hand side of Eq. H23|l is simply proportional to 
the intrinsic entropy perturbation of the scalar field F, 
hence in the large-scale limit this is forced to vanish if <i> 
is to remain small. This is a known result already shown 
in Refs. 

By contrast, once a fluid is added we have 



CgSptot is the non-adiabatic pres- 



where Jpnad = <5ptot 
sure perturbation. Note that even on large scales TZ 
can evolve due to the presence of a non- vanishing non- 
adiabatic pressure perturbation, as recently stressed in 
different works [Tsl . The non-adiabatic pressure per- 
turbation depends on the intrinsic and relative entropy 
perturbations |l6l [TtI I , and in our case we find 



Sp 



'nad 



Ptot 



[{wt - c^,^)S + {1 - cl^)T] 



(27) 



From now on we find it convenient to describe the system 
in terms of the gauge-invariant variables <&, TZ, S and F, 
rather than the set of variables <&, S[, S^, and (5^. Note 
that such a change of variables is completely determined 
by Eqs. H2U|) . H21|l and the expression 



37 



^5^ + rti5{) — 



2 



(28) 



obtained from Eqs. ^ and (|25|l . In the next section 
we find a first-order system of dynamical equations ex- 
pressed in these variables. 



2^2 



[n^T + ili {Si - 3i/7fVf )] 



(24) 



This equation shows that the fluid comoving density per- 
turbation can compensate the scalar field intrinsic en- 
tropy perturbation, and, as a result, in the presence of 
a fiuid it is possible to have a non-vanishing scalar field 
intrinsic entropy perturbation even on large scales. Note 
that the presence of the fluid changes the structure of 
the equations even if it is a sub-dominant component of 
the total energy density. This is because the fluid creates 
a new set of hypersurfaces, those on which its density 
is uniform, which need not align with hypersurfaces of 
uniform scalar fleld density. 

Since we are interested in studying the evolution of 
isocurvature and adiabatic modes, we flnd useful to use 
the gauge-invariant comoving curvature perturbation 



_ 2(iJ$ + $) ^ 
7^= ^ 

The equation of motion for TZ is given by [ll 



TZ 



— $ + 47rG(5pnad 



(25) 



(26) 



III. MATRIX FORMULATION 

A. Evolution equation 

In the long- wavelength limit (k/aH ^ 1) Malquarti 
and Liddle 4_1| were able to express the dynamical equa- 
tions in a flrst-order matrix formulation, using N = 
log(a/ao) as a time variable. They took as basic vari- 
ables $, 5i, d^, and (5^. Here we show that our set of 
variables can bring the matrix into an even more effi- 
cient form. Moreover, we compute the general equations 
without the long-wavelength approximation. 

We define the vector v = (<i>, TZ, S, F)^ and use a prime 
to denote a derivative with respect to A''. Lengthy but 
straightforward algebra leads to the expression 



Mo + Mi 



'A 1 



X V 



(29) 



where the only relevant matrix for the long-wavelength 
approximation [k / aH = 0) is given by 



/-(I + 37/2) 37/2 

n^iwf 

3(w^ 

\ -37/2 



)h 



)/7 



Wi) + i-ii^i{wi - ci^)h 37frif (1 - c2^)/7 

3(u;^ - 7/2) ) 



(30) 
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and the two matrices incorporating spatial gradients are 



Ml 



and 
























" ] 


V 








1/3 


1/3 







-1/3 


-1/3/ 



M2 








27^/97 
-27^/97 0. 



(31) 



(32) 



S' ~ and therefore on large scales F decays exponen- 
tially^ with decay rate —3 + 37^/2, dynamically recov- 
ering the single scalar field case for which the intrinsic 
entropy perturbation vanishes (cf. Section III CI) . 

When the fluid is completely absent, the matrices in 
Eq. (|29() reduce to 3 x 3 matrices, in the variables $, TZ 
and r, but the constraint in Eq. (|23|l allows to eliminate 
one more degree of freedom. For example, using Eq. I|23|) 
one can find 2x2 matrices for TZ and F, or if one addition- 
ally goes to the large-scale limit, the constraint equation 
Eq. (|23|l forces F to vanish and gives A^o as a 2 x 2 matrix 
for $ and TZ. 



The different non-vanishing entries clearly show the 
couplings between adiabatic and relative/intrinsic en- 
tropy perturbations on large scales (Ado) and on small 
scales {Ml and AI2). The first two fines of the matrices 
(dynamical equations for <I> and TZ) are straightforward 
from Eqs. H26I) and (|25() . The equation for the relative 
entropy, here expressed as 



Adiabatic condition 



S' 



3{wtp — Wf) + 
37ff^t(l-c2 ), 



37f»f (lUf - cg^) 



(33) 



7 



3 3 



27V 
9 7 



has been obtained both in the context of multiple inter- 
acting fluids 16, 17, 2QJ, and in the framework of infla- 
tion when several interacting scalar fields are present [3 j 
Eq. (|33|l being a particular case. However, in general it is 
not possible to find a dynamical equation for the intrin- 
sic entropy perturbation of a given component without 
knowing its underlying physics. In the case under study 
we are able to fully specify the evolution of the system 
through the equation for the intrinsic entropy perturba- 
tion of the scalar field as 



+ 



-l^TZ 

27^ 
9 7 



7 

1 



-S- 



(34) 



Eqs. 1)33(1 and ()34|l show that, on large scales, the relative 
entropy perturbation and the intrinsic entropy perturba- 
tion of the scalar field are mutually sourced and evolve 
independently of the curvature perturbations. In partic- 
ular, Eq. ((34(1 confirms the conclusions drawn from the 
constraint equation Eq. 1(24(1 . namely that in the presence 
of the fluid it is possible to have an intrinsic entropy per- 
turbation relative to the scalar field even on large scales. 
When the fluid is very sub-dominant (fif ~ 0), we have 



The adiabatic condition requires that the relative en- 
tropy perturbation S and the intrinsic entropy pertur- 
bation F vanish. From our equations it is immediately 
clear that on large scales {k/aH ^ 1, so that only Mo 
need be considered) if the perturbations are initially adi- 
abatic then S and F remain zero for all times. In this 
case TZ is constant and $ rapidly approaches its asymp- 
totic value $ = i'jTZ/ (2-1-37) (for constant or sufficiently 
slowly varying 7). This demonstrates that adiabaticity 
on large scales holds regardless of any time-dependence 
of the background variables such as and c^^ (this was 
already pointed out in Ref. ^3). In fact, preservation of 
adiabaticity is implied by the separate universe approach 
to large-scale perturbations jl5|. However adiabaticity 
will be broken once the perturbations move out of the 
large-scale regime, with the matrices Mi and M2 sourc- 
ing S and F through the curvature perturbations TZ and 
$. In particular, note that this is also true for the single 
scalar field case, as is evident from looking at Eq. 1(23(1 . 
On the other hand, as we will see, if an isocurvature per- 
turbation is initially present it can be wiped out on large 
scales by the scalar field dynamics. 

Aspects of these results have appeared in previous 
works 0, 0) ^> ^1 J but without noting that adiabaticity 
is always preserved on large scales. Our set of variables 
makes unambiguously clear the fact that the adiabatic 
condition is not an instantaneous notion on large scales, 
and holds independently of the evolution of the back- 
ground. 



^ Actually, r remains constant in the special case = 1, but 
when the field is dominant this equation of state is usually not 
considered and anyway would rapidly evolve towards < 1. 
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c2 


Kinetic 





2 


1 


Potential I 








1 


Potential II 








-2 - wt 


Usual Tracker 





If 




Perfect Tracker 




7f 





TABLE I: Values of the three parameters 57^, 7^, and c^^ 
during the four different possible regimes of a quintessence 
scenario (a particular scenario would feature only one type of 
tracker regime) until the field starts dominating. 



IV. APPLICATION TO THE QUINTESSENCE 
SCENARIO 

A. Analytical description 

In this section, we discuss the large-scale evolution of 
perturbations in quintessence scenarios. As described in 
the Appendix, before the quintessence field starts dom- 
inating the evolution of the universe its dynamics can 
feature up to four different regimes during which the co- 
efficients of the matrix A^o are constant. These are sum- 
marized in Table m Note that, as compared to Refs. [sj 
and pd| . we altered the names of two regimes (poten- 
tial I and II) in order to make our explanations clearer. 
Now, following Ref. for each regime it is possible 

to perform an eigenvector decomposition of the matrix 
A4o and therefore compute analytically the large-scale 
evolution of the perturbations (i.e. during each regime 
V can be written as a sum of four terms proportional 
to Viexp{niN) = Via"' for i = 1 to 4 respectively). 
However, the matching conditions between the different 
regimes are not obvious as S and F contain non-trivial 
functions of the background. In that respect, the formu- 
lation in Rcf. 11] is more appropriate when following the 
modes over different regimes, since to a first approxima- 
tion $, 6f, 5ip, and can be taken as conserved through 
the transitions between regimes. 

First of all, it is easy to find that for any regime M.q 
possesses two eigenvectors 



VI = (37,2 + 37,0,0) 
V2- (1,0,0,0) 



ni = 0, 

n2 =-1-37/2, 



(35) 



where rix is the eigenvalue of Vx . These two vectors cor- 
respond to the two well-known adiabatic modes, the first 
one being constant, and the second one rapidly decaying. 

Now, it is possible to find the two remaining entropy 
modes for a general case, but they cannot be expressed 
in a simple form. Nevertheless, it is straightforward and 
more clear to perform an eigenmode decomposition by 
considering each regime separately. The modes are given 
for each regime in Table ^1 For simplicity, we do not 
display complicated expressions; for a detailed analysis 
the reader should refer to Ref. [ll| ■ Most of these results 
have already been discussed in that paper, but here we 



Kinetic regime 
V3 = (0,0,0,1) 



ns = 3 — 37f/2 



V4 = (0,0,2 - 7f,7f) 714 = 



Potential regime I 

V3 = (0,0,0, 1) n3 = -3-37f/2 

V4 = (0,0,2 - 7f,7f) n4 = -6 



Potential regime II 

V3 = (0,0,-2 -7f,7f) n3 = 

V4 = (0,0,-2, 1) n4 = - 



3 + 37f/2 



Tracker regime (7^ < 7f) 

V3 = . . . Re(n3) < 

V4 = . . . Re(n4) < 



TABLE II: Eigenvectors and corresponding eigenvalues of the 
matrix A4o of Eq. 1291 according to the different regimes. 



would like to comment further in the light of our new set 
of variables and new findings. 

As is well known, we see that entropy perturbations 
decay during the tracker regime ,3, illj . This is due 
to the scaling and attractor properties of that regime. 
More striking is the fast-growing mode during the kinetic 
regime and the constant mode during the second poten- 
tial regime. This contradicts the claim by Brax et al. Q 
that the final value of the quintessence perturbations is 
insensitive to the initial conditions. The difference comes 
from the fact that these authors considered the basic vari- 
able dip and its time derivative. They found that there are 
two decaying modes for every regime. However, this does 
not mean that observationally-relevant variables (such as 
6^p) are decaying, since one has to take into account the 
evolution of background quantities (such as p^) as well. 
Our set of variables is therefore more appropriate. 

Since in general we expect that there could be an ini- 
tial relative entropy perturbation S (for example in the 
case of a quintessence field present during infiation j23j ) 
and since 5" sources F, we can expect a non-zero intrinsic 
entropy on large scales which would then evolve accord- 
ing to our set of equations Eqs. and Now, using 
the results given in the Appendix, Eqs. (p??^ and (|^ . we 
can show that the growth of F during the kinetic regime 
— exp(AAfk?^3(k)) = exp(37t AAfp/2) — is exactly com- 
pensated by its subsequent decay during the potential 
regime I — exp(AAfpin3(pi)) = exp(— 37t AAfp/2). As 
a result, after the three regimes preceding the tracker 
regime, entropy modes are neither enhanced nor sup- 
pressed. However, as shown in the Appendix, accord- 
ing to the initial conditions, the potential energy of the 
field can undergo a very large drop before it reaches a 
constant value and the kinetic regime starts. In general. 
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this transition phase could last a non-negligible number 
of e-foldings and would feature the same eigenmodes as 
the kinetic regime, in particular the same growing mode. 
Since this first phase of growth would not be compen- 
sated by the decay during the first potential regime, there 
remains the possibility that entropy modes may be en- 
hanced at the beginning of the tracker regime. Now, 
if the tracker regime is not long enough to erase com- 
pletely these entropy perturbations by the time the field 
becomes non-negligible, we may be able to see an imprint 
of these initial perturbations in observations As a 

result, we can see that in a quintessence scenario the ini- 
tial conditions, as well as the history of the evolution of 
the background, are relevant when considering the late- 
time value of the perturbations. 

In this respect many different assumptions can be 
made. In Ref. 21], Kneller and Strigari assumed equipar- 
tition as an initial value for the field, and in most cases 
this led to a field dynamics featuring a very long kinetic 
regime followed b y th e potential regimes and no tracker 
regime. In Ref. |23|, de la Macorra studied an actual 
physical model in which the quintessence field is a dark 
condensate which arises after a phase transition. Its evo- 
lution starts in the kinetic regime and again does not 
feature a tracker regime. In both scenarios entropy per- 
turbations would still be present today. Alternatively, 
Malquarti and Liddle '23'| studied the evolution of a 
quintessence model during inflation in order to investi- 
gate the initial conditions of the quintessence field at the 
beginning of the radiation-dominated era. They found 
that typically the tracker starts at low redshift after a 
long period of potential regime II, but again, as a result, 
entropy perturbations generated during inflation could 
still lead to observable consequences. Finally, note that 
in general the initial entropy perturbations do not need 
to be of the order $ and may be much larger. 



B. Numerical examples 

In order to illustrate our results, we carried out simu- 
lations for two quintessence models in a realistic universe 
(with radiation, dark matter and dark energy). 

The first example is an inverse power-law model 3j 
V{(p) = Voi(p/mpi)-°', with a = 1 and = IQ-^^^ m^.^, 
starting with the initial conditions at = — 50 given by 
^ 10~^"pf, T = V and S = -T = E-,. The results are 
shown in Fig. ^ where we display the evolution of some 
background variables and of the two entropy variables S 
and F. After a very sharp transition towards the kinetic 
regime, the field undergoes the four regimes described in 
TableQ]before starting dominating. First, we can clearly 
identify these regimes and see the growing, constant and 
decaying behavior of S and F corresponding to the modes 




-20 I ' • ' • • • ' 1 

-50 -40 -30 -20 -10 

N 



FIG. 1; Evolution in a realistic universe of background quan- 
tities (top) and entropy perturbation variables (bottom) of 
a quintessence field undergoing the four different possible 
regimes before its domination. We use an inverse power-law 
potential V{ifi) — Vo{}p/m-p\)~'^ (see Section FlV BI for parame- 
ters). Note the transition between radiation domination and 
matter domination at A'^ ~ —9. 



displayed in Table HTl^ We note the oscillations during 
the tracker regime and the non-trivial evolution through 
the transitions (yet keeping the same order of magni- 
tude). We also observe that in this particular case the 
transition phase towards the kinetic regime is too short 
to enhance the entropy mode significantly, leading to the 
result that at the beginning of the tracker regime S and 
F have about the same amplitude as at the initial stage. 
In other words they are neither enhanced nor suppressed 
by the dynamics of the field, until tracking begins. 

The second example is a double-exponential model |3| 
V{Q) Vb[exp(-Q;K(^) -I- exp(-/3K(/?)], with a = 1000, 
(3 — 1, Vq = 10~^^^TOp[ and k — y^Stt/Stopj, starting 
with the initial conditions at A^ = — 50 given by p^p = 
5 X 10-3pf, T = y and 5 = and F = 0. The results 
are shown in Fig. |2 We display the same variables as 
in Fig. ^ In addition, in order to observe the transition 
phase preceding the kinetic regime, we also display V' /V ^ 
where V is the potential energy density of the field — as 
described in the Appendix, the kinetic regime starts when 
drops under '-^ 1. First, we note that although 
F initially, it evolves very quickly to be of the same 
order as 5'; as previously explained, this is due to the 
coupling between S and F. Then, we observe the same 
behavior as for the first example, but this time we can 
see that the transition phase towards the kinetic regime 
lasts a few e-foldings. As a result, the entropy modes at 



^ Note that contrary to what the modes displayed in TablelfTlsug- 
gest, S is constant during the potential regime I. This is because 
r is many orders of magnitude larger than 5, and 1 — Cg^ is not 
exactly 0. 
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N 

FIG. 2: As Fig.0 but now for a double-exponential potential 
V{Q) = Vb[exp(— QK<^) + exp(— /3k<^)] (see Section liV Bl for 
parameters). We also display V' /V (to be read on the right- 
hand side of the graph). 



the beginning of the tracker regime are nearly a factor 
100 larger than initially. This shows that, as discussed 
in Section llV Al entropy modes can be enhanced before 
the field reaches the tracker regime. 

We have shown that, in general, one must take into 
account the initial conditions for the quintessence field 
and its perturbations in order to make any prediction. 



V. DISCUSSION 

We have explored the nature of scalar perturbations 
for a universe filled with both a scalar field and a per- 
fect fluid. We have introduced a useful set of variables 
and have provided a full analysis including spatial gradi- 
ents. In particular, we have focused on the isocurvature 
perturbation modes and on the degrees of freedom which 
completely characterize the system. While for the case 
where only the scalar field is present its intrinsic entropy 
perturbation is forced to vanish at linear order for su- 
perhorizon scales, the presence of a fluid — even if sub- 
dominant — allows the possibility for such an intrinsic 
contribution to be present on large scales. However, in 
the case of a very sub-dominant fluid the intrinsic en- 
tropy of the scalar fleld decays, dynamically recovering 
the single scalar field situation. 

We have recast the basic evolution equations in a 
rather simple matrix formalism in terms of the gauge- 
invariant variables for the adiabatic and isocurvature per- 
turbations, taking into account the dynamics of the per- 
turbations when a given wavelength re-enters the hori- 
zon. In particular, we have obtained an equation for the 
intrinsic entropy perturbation which shows that, on large 
scales, an initial adiabatic condition is indeed preserved, 
regardless of the evolution. Only when the perturbations 
approach the horizon are the adiabatic and entropy per- 



turbations fully coupled together. 

Finally, we have applied our formalism to the 
quintessence scenario. In this case we have analyzed the 
large-scale evolution of the adiabatic and entropy pertur- 
bations in the different regimes which the quintessence 
scalar field dynamics may feature. As is well known, 
entropy perturbations are suppressed durin g th e track- 
ing regime, but, as already shown in Ref. [T^l, during 
the kinetic regime one entropy mode undergoes an expo- 
nential growth. We have shown that it is then exactly 
compensated by an exponential decay during the first po- 
tential regime and then remains constant during the sec- 
ond potential regime. Therefore, after the three regimes 
preceding the tracker regime, entropy modes are neither 
enhanced nor suppressed. However, we discussed the re- 
maining possibility of an enhancement during the short 
transition phase preceding the kinetic regime. We have 
studied two special cases numerically and have confirmed 
our analytical analysis. Moreover, we have observed that, 
in one of the cases, at the beginning of the tracker regime, 
entropy perturbations are larger than initially, and there- 
fore we have concluded that entropy mode enhancement 
is possible. 

To summarize, we have shown that in general it is 
incorrect to assume that the observational imprint of 
quintessence perturbations will be independent of the ini- 
tial conditions, because entropy perturbations can still 
be present when the quintessence energy density is no 
longer negligible. Note that this can happen as soon as 
the tracker starts, long before quintessence domination. 
In this case, entropy perturbations would feed curvature 
perturbations, but then would slowly decay to become 
negligible today. 
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APPENDIX: DYNAMICAL REGIMES OF A 
QUINTESSENCE FIELD 

As described in Ref. 0, a tracking quintessence field 
can feature up to four different dynamical regimes when 
in presence of a dominant fluid with constant equation of 
state. Here we clearly demonstrate the existence of these 
regimes and compute some relevant parameters. We will 
use the same notation and definitions as in the main body 
of the article. 
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FIG. 3; Example of a quintessence field evolution featuring 
the four possible regimes before domination: kinetic, potential 
I, potential II and tracker. Note the transition phase towards 
the kinetic regime. The different parameters are described in 
the Appendix. 



We assume that the quintessence field features a track- 
ing solution pti{N) — note that it does not need to have a 
constant equation of state. This means that at each time 
(i.e. for each value of H) there exists a stable field con- 
figuration for which its kinetic energy density T = <yj^/2 
and its potential energy density V are of the same or- 
der, and hence in Eq. the 'friction' term due to the 
Hubble expansion and the slope of the potential balance 
each other. We will show that, according to the initial 
conditions, the scalar field can feature up to three dif- 
ferent regimes before it reaches the tracker. We assume 
that the field is always subdominant and therefore does 
not influence the evolution of the universe, especially the 
evolution of H. As a result we have H cx exp(— 37f7V/2). 
Using Eq. Q it is easy to see that 



Y 
YL 



dV/Aip ypXT 



V 



H 



(36) 
(37) 



In order to help the reader to follow the explanation, 
an example of a quintessence field evolution featuring 
the four possible regimes before domination is shown in 
Fig.E 

Let us start by looking at the initial condition T < 
V » ptr at time TVj. The field, and hence the slope 
of the potential, has the same value as for the tracker 
configuration at an earlier time < Ni, but because 
H{Ni) ^ i?(iV_) the friction term is actually negligible 



and the field fast-rolls down its potential {V <C T') and 
its kinetic energy almost instantaneously dominates its 
energy density. At some time TVk shortly after iVj (at the 
latest when T Ptr), the potential freezes at some value 
V+ <C Ptr corresponding to the tracker configuration at 
a later time N+ > TVk, and, since H{N) > H{N+) for 
N < 7V+, it remains frozen until N^. In this case, the 
evolution of T can easily be computed analytically. We 
assume that at the tracker solution has an equation 
of state 7t"^, express dV/dip as a function of V+ and 7^"^, 
and for A^k < ^ < ^+ we find 



T(m- - ltr)V+ ^, 



(38) 



where C is a constant of integration depending on the 
initial conditions. 

Let us explain this behavior. Starting from the time 
A'k the field is in the kinetic regime: = T (x. exp(— 67V) 
and V = constant. At some time A^pi the field reaches 
the configuration T ~ y, and since the friction term is 
still extremely large {H{Npi) H{N+)) the kinetic term 
keeps on decaying and the field enters the first potential 
regime: p^ — V — constant and T oc exp{—6N). At 
some time iVpn the term V' /T can balance the friction 
term and T starts growing again. Here begins the second 
potential regime: Pip = V = const and T cx exp(37fiV). 
Finally, at iV_|. the field enters the tracker regime: T ~ 
V ^ Ptr- Note that the whole evolution described above 
goes through all the possible initial conditions. 

Now we can compute a few parameters. First, using 
the solution for T given in Eq. (|38|l and noting that c^^ = 
— 1 — T' /3T, it is straightforward to recover the values 
displayed in TableU In addition, let us define AN^ as the 
number of e-foldings that the field spends in the regime 
'r' and also AA^p ee AA^pi -I- AiVpn. Using Eq. ||2Zll, the 
fact that at iVk and N+ we have V /V = -37^+ -1 
and the evolution for T and H described above, we find 



AA^k 



7f 



If 



■AiVp, 



7f 



7f 



ANp, = ^AA^pii^ 

z ^ + 7f 



AN„ 



(39) 
(40) 



These last two results are used in Section FlV Al to show 
that the growth of one of the entropy modes during the 
kinetic regime is exactly compensated by its subsequent 
decay during the first potential regime. 
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